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On Transient Solutions of the "Baffled Piston" 

Problem 

F. Oberhettinger l 

(September 14, 1960) 

The acoustic field produced by the movement of a piston in an infinite rigid wall for 
arbitrary time dependency of the motion is given. 

1. Introduction 

The case of the time harmonic movement of a piston membrane in an infinite rigid wall 
("baffle") can easily be generalized to the case when the motion of the piston is not periodic 
but an arbitrary function of time. Such transient solutions have become of considerable in- 
terest in recent times (for a detailed treatment of the propagation of such sound pulses, see [5] 2 ). 
The procedure for the case treated here is the same as used elsewhere [7], i.e., the Green's func- 
tion for the exponential decay case (modified wave equation AU— y 2 U=0, y=ik) is used to 
obtain the solution for the pulse problem. The accustic field (velocity potential) for the lime 
harmonic movement of the piston include representations given by Bouwkamp [1], King [6], 
and Wells and Leitner [9]. The first of these contributions gives the solution in the form of a 
series expression while the second and third involve integral representations that are obtained 
using integral transform methods (Hankel transform [4, p. 7i]] and kehedev transform [4, p. 75] 
respectively). These representations can he used to treat the general case of an arbitrary move- 
ment of the piston. In view of the method to be employed here, such representations should be 
used for which the inverse Laplace transform of the velocity potential with respect to the purely 
imaginary wave parameter y= ik can be given. Such an expression can be obtained in a direct 
way by regarding each point of the moving disk as an accustic point source and integrating 
over all points of the disk. 

2. Exponential Decay and Time Harmonic Solution 

The piston is represented by an infinitely thin circular disk of radius a located in the x,y- 
plane with its center at the origin of a three-dimensional Cartesian system of coordinates. The 
remaining part of the x,y-jAane consists of a rigid wall. The movement of the disk is time 
harmonic of the form 

v z =e tut (1) 

where v z is the velocity of the disk. 

Each point Q of the disk can be considered as a point source with the velocity potential 

$=-,^«'"VW, v=PQ, (2) 

(a)=fcc, c is the velocity of sound), which represents Green's function of the wave equation in 
free space (see fig. 1). 
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Figure 1. 



Upon substituting 



k= — iy, cc=—iyc, 



(3) 



one obtains, after isolation of the time-dependent factor e yct for the free-space Green's func- 
tion of 

Au-y 2 u=0, (4) 

the expression 



u =-y4^p'd<p' d p f , 



r=[p 2 +p> 2 -2pp f cos (^-V) + 2 2 ] 1/2 . 



(5) 



(6) 



The transition from k to y by (3) amounts to the transition from a wave problem to an expo- 
nential decay problem. At first the exponential decay problem for positive real y is solved and 
the solution for the wave problem can be obtained by returning from y to k. This, however, 
is not necessary here, since the exponential decay solution will be used in order to solve the 
problem of an arbitrary moving disk. If the expression in (5) is integrated over the whole disk, 
one obtains — observing (6) : 

/*2ir f*a 

C7= -f [p 2 + P ' 2 -2pp'cos(^') + 2 2 ]- 1/2 

• exp {[p 2 + P > 2 -2pp' cos 0-V) + 2 2 ] 1/2 } p'dp'dv'. (7) 

Here t\ is the "density distribution" over the disk. Since the disk moves as a whole, this prop- 
erty is a constant; it is determined by the fact that according to (1) 



UJ 2=0 io! 



for p<Ca, 
for p^>a. 



(8) 



This function U as defined in (7) satisfies 

(a) AU-y 2 U=0, 



(b) 
(c) 






■■— 1 for p<a, 



W_ 
dz 



=0 for p~>a at 2=0, 



lim U-^Q when p or |z|-><». 



In order to perforin the integration in (7), one replaces the integrand in (7) using the well- 
known formula [3, vol. 1, p. 191, 9], 

pJo[X(< 2 -7 2 ) 1/2 ]e- z '^=(3 2 +X 2 )- 1/2 e-^ z2 + x2)1/2 , 2^0 

or, upon substituting 

^- 7 2= T 2 j t 2 =T+y 2 

(2 2 +X 2 )- 1/2 e- 7(z2+x2)1/2 = r" t / (XT)T(r 2 +7 2 )- 1 ' 2 e- z( ' 2 + 1 ' 2)1/2 rfr (9) 

If this formula is used to replace the integrand in (7), one obtains upon interchanging the 
order of integration 

Z7=-^J o %(T 2 + 7 2 )- 1 / 2 e-^ 2 +^ 1/2 x/J o a J o 2 Vo(r[p 2 +p' 2 -2p/cos(^-/)] 1 / 2 )p'rfp'^'"ydr. 

(10) 

To evaluate the integral in (10), one uses the addition theorem for the Bessel functions [4, 
p. 101] 

J n {T[p 2 +p /2 -2pp' cos (ip— <p')\ 1/2 } = ^>2 e m J m (Tp)J m (T P ') cos[m(<p-<p')\, 

e =l, e m =2, m=l,2,3, . . . 

Integrating with respect to <p' and observing that only the term m=0 gives a contribution, 
yields 

c7=-^J"r(r 2 + 7 2 )-" 2 e -^ 2 +^ 1/2 J (rp)|j'yj l ,(rp')'/p , }rfr. 
Finally, using the result [4, p. 45, 1] 

f «7„(Tp')p'<*p'=-«7i(Ta), 

J() T 

one obtains 

^=-^r(r 2 +70- 1/2 ^ 2(T2+72)1/ Vo(rp)J 1 (ra)dr. (11) 

Now, according to (1) 

-(SLr 1 ' p<a > -(SLo =0 ' p>a - 

therefore by [3, vol. 2, p. 14], 

VcW 2 =o 2 ar? \0, P >a. 
Hence 

V=-2. 

The final representation becomes therefore 

U = a f"( T2 +y 2 )~ l/2e ~ 2(T2+y2)m J o(rp) Ji {ra)dT, (12) 

Clearly U satisfies the conditions a, b, c, stated following (8). The restriction z^0 in (12) 
means no loss of generality, since it is sufficient to consider the upper half space. If one wishes 



to return to the wave problem one has to replace y by ik according to (3) and one arrives in 
this case at King's solution [6]. 

Approximate evaluations of King's and related integrals for various cases have been 
carried out [2, 6, 8]. 

3. Transient Solutions 

The expression (12) is now used to investigate the corresponding pulse problem. Let the 
movement of the piston be represented by a pulse starting at £=0; let the pulse function be 
g(t), and g(t) = Q for £<0. Under the assumption that the pulse function g(t) is such that it ad- 
mits its representation by Laplace's integral formula 

g{t) = {2iri)-' ( C [ P g^e^dr'] e^dy, (13) 

one can represent the expression for the acoustic field generated by the pulse in the form 

^(t) = (27ri)- 1 ( C [{ g(T)e- y TdT~\ue yt dy (14) 

Jc-ioo LJo J 

where U is the representation (12) for the " exponential decay" field. (The velocity of sound 
has been chosen to be unity.) 

If the order of integration can be interchanged one obtains 

$(t)= r g(r)[(2iri)- 1 f ^ Ue nt - T) dyldr. (15) 

•7 L J C — lca J 

The integral in the bracket is the inverse Laplace transform of the " exponential decay" solution 
with respect to the wave parameter y. If this is known, then the velocity potential <£(£) of the 
acoustic field for the case of an arbitrary movement of the piston can be obtained by a further 
integration involving the pulse function g(t) according to (15). The inverse Laplace trans- 
form of the integrand in (12) with regard to y is known [3, vol. 1, p. 248]: 

(art)- l T °° (r 2 +7 2 )" i 6-'^ 2 +r 2 )^ e r^ T= J r Jo[r(^-S 2 )*] (16) 

Jc-ioo V* ^ 

according as 0>z or t<^z, respectively. 

The acoustic velocity potential due to a pulse function g(t) ("Dirac" pulse movement of 
the piston at t=0) is by (14) 

^ D (t) = 2iri)- l f C l "u o yt dy. (17) 

J C — ioo 

Hence, (15) can be written as 

/•CO 

•(0—1 g(r)*D(t-T)dT. (18) 

Therefore, by (12) and (16) 

' - l {ra)J (rp)J ii [r{t'-Z^\dr, t>Z. (19) 



'»«>={«/> 



The integral is known. One has [3, vol. 2, p. 21]: 

r 0, 0<y<\b-c\ 

Jo(bx)J (cx)J l (xy)dx=< (iry) ~ l arc cos f - , y b \b—c\<Cy<b+c 

<y~\ y>b+c 



(20) 



Iii order to evaluate (19) identify y with a, b with (P—z 2 ) 1 ' 2 and c with p. Two different cases 
have to bo distinguished: 

Case A. a>p, 

Case B. a<p, 

i.e., the two possible different cases where the projection S ot the point of observation P into 
the £,?/-plane lies either inside or outside the disk are treated separately (fig. 2 and fig. 3 
respectively). 
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Figure 2. 



Figure 3. 



The largest and the smallest distances of a point on the circumference of the disk from the 
point of observation P are denoted by R and R f respectively. 



R=[z 2 +(a+p) 2 ] 1/2 , R' = [z 2 +(a-p) 2 Y /2 . 

One obtains then for $ D (f) by (19) under consideration of (20) and (21): 
For the case A, a>p: 



(21) 



**(*) -o, 



t<z 



For the case B, a<Cp: 



=1, z<t<R' 

r^2_ g 2_j_ 2_ a 2-| 

=x -1 arccos 1 2 (? 2 — g 2 )'/ 2 J' #'<K-K 

=0, *>#. 

*j)(*)=o, KB' 

-i fi 2 — 2 2 +p 2 — a 2 "| p//,/p 

=tt arccos [_~2 p ( ^_ g2) 1/2 j' « <K« 

= 0, <>fi. 



(22) 



(23) 



The physical interpretation of these results is obvious. The parts in (22) and (23) that 
are different from zero represent the transient acoustic field at a point P due to a "Dirac" 
pulse movement of the piston for a time which lies between the smallest and the largest distance 
of the point of observation P from the disk. Outside of this the field is zero. (Note that the 
velocity of sound was assumed to be unity.) The "Dirac" pulse solutions (22) and (23) can 
now be used to construct the acoustic field in case of an arbitrary time dependency g(t) of the 
movement of the piston according to (18). One obtains immediately: 

For the case A, a>p: 

m= ^n' g{T) arcc ° s { %[t-ry-%^ } dT+ fz 9 {r)dr - m 

For the case B, a<p: 

•®=*"0 (T) arCC ° S { %ll-r^-%- } ^ (25) 

The case p=0 (point of observation on the axis of the disk) reduces to 

•O-Ji^.^W*. (26) 

It has to be remembered that g(t)=0 for t<C0. Therefore no possibly negative part of the 
interval of integration in (24) to (26) gives a contribution. Furthermore, the time parameter t 
in the above formulas has to be replaced by ct, where c is the velocity of sound. 
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